Stabilizing a purely dipolar quantum gas against collapse 
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We report on the experimental observation of the dipolar collapse of a quantum gas which sets in 
when we reduce the contact interaction below some critical value using a Feshbach resonance. Due 
to the anisotropy of the dipole-dipole interaction, the stability of a dipolar Bose-Einstein condensate 
depends not only on the strength of the contact interaction, but also on the trapping geometry. We 
investigate the stability diagram and find good agreement with a universal stability threshold arising 
from a simple theoretical model. Using a pancake-shaped trap with the dipoles oriented along the 
short axis of the trap, we are able to tune the scattering length to zero, stabilizing a purely dipolar 
quantum gas. 



Interactions between atoms dominate most of the prop- 
erties of quantum degenerate gases p^l . In the ultracold 
regime these interactions are usually well described by 
an effective isotropic zero-range potential. The strength 
and sign of this contact interaction is determined by a 
single parameter, the scattering length a. The contact 
interaction is responsible for a variety of striking prop- 
erties of quantum gases. Strongly influencing the excita- 
tion spectrum of the condensate it gives rise to e.g. the 
superfluidity of Bose-Einstein condensates (BEG) or the 
existence of vortex lattices. The contact interaction also 
plays a crucial role in the physics of strongly correlated 
systems like in the BEG-BGS crossover Q or in quantum 
phase transitions like the Mott insulator transition [s!]. 

Another fundamental topic is the question of the exis- 
tence of a stable ground state depending on the modulus 
and sign of the contact interaction. In the homogeneous 
case repulsive contact interaction (a > 0) is necessary for 
the stability of the BEG. In contrast, if the contact inter- 
action is attractive (a < 0), the BEG is unstable. This 
instability can be prevented by an external trapping po- 
tential. The tendency to shrink towards the center of the 
trap is in that case counteracted by the repulsive quan- 
tum pressure arising from the Heiscnberg uncertainty re- 
lation. Detailed analysis 0] yields that a condensate is 
stable as long as the number of atoms N in the conden- 
sate stays below a critical value A'crit given by 



(1) 



where aho is the harmonic oscillator length and fc is a 
constant on the order of 1/2. This scaling, as well as 
the collapse dynamics for N > Ncrit, have been studied 
experimentally with condensates of ^Li and ^^Rb 0, 
i- In[i,[l3 

the atom number dependance of the collapse 
of mixtures of bosonic ^^Rb and fcrmionic ^"^K quantum 
gases has been investigated. 

Being anisotropic and long-range, the dipole-dipole in- 
teraction (DDI) differs fundamentally from the contact 
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interaction. Besides many other properties, the stabil- 
ity condition therefore changes in a system with a DDI 
present. Gonsidering the case of a purely dipolar con- 
densate with homogeneous density polarized by an ex- 
ternal field, one finds that due to the anisotropy of the 
DDI, the BEG is unstable, independent of how small the 
dipole moment is 11|. As in the pure contact case a 
trap helps to stabilize the system. In the dipolar case, 
however, it is not only the quantum pressure that pre- 
vents the collapse but more importantly the anisotropic 
density distribution imprinted by the trap. Consider a 
cylindrically symmetric harmonic trap 

Vtrap{r,z) = im (w^r^ 



(2) 



with the dipoles oriented along z and r being the distance 
from the symmetry axis. As can be seen in Fig. lA, in 
a pancake-shaped trap (aspect ratio A — ujzl^r > 1) the 
dipoles predominantly repel each other and the BEG is 
stable. In contrast, a cigar-shaped trap (A < 1, Fig. IB) 
leads to mainly attractive forces and hence to a dipo- 
lar collapse. Following this simple argument one expects 
that in the prolate case a positive scattering length a is 
needed to stabilize the BEG, whilst in the oblate case one 
can even afford a slightly negative a. The dependance of 
the stability of a dipolar BEG on the trap aspect ratio A 
and scatte ring len gth a has been extensively studied the- 
14| , and is experimentally investigated 



oretically l^, 13, 
in this paper. 

Our measurements are performed with a BEG of 
^^Gr [iBj which is to date the only experimentally ac- 
cessible quantum gas with observable dipole-dipole in- 
teraction [3, 113] ■ To compare contact and dipolar in- 
teractions we introduce the length scale of the magnetic 
DDI 



Odd 



(3) 



The numerical prefactors in add are chosen such that a 
homogeneous condensate becomes unstable to local den- 
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FIG. 2: The atom number A'' in the condensate as a function 
of a for two traps having different aspect ratios A. The solid 
lines are fits to Eq.[5]used to determine the critical scattering 
length acrit (see text). 



FIG. 1: (A, B) Intuitive picture of the trap geometry de- 
pendance of the BEG stability: In an oblate trap the dipoles 
mainly repel each other, whereas in a prolate trap the interac- 
tion is predominantly attractive. (G) The different trapping 
geometries are realized by the crossed optical dipole trap (red) 
and an additional ID optical lattice (green). The magnetic 
field is pointing along the symmetry axis z of our traps. 



sity perturbations for a < add uM- Chromium has a 
magnetic dipole moment oi fi = 6fj,B (fJ-B the Bohr mag- 
neton), add — 15ao, where ao is the Bohr radius. Far 
from Feshbach resonances, a takes its background value 
abg — lOOao 19] and the BEC is stable for any A. To 
explore the unstable regime we thus reduce the scatter- 
ing length, which in the vicinity of a Feshbach resonance 
scales like 



abg 



AB 

B-Bn 



(4) 



with the apphed magnetic field B. To be able to tune 
a accurately we use the broadest of the resonances in 
^^Cr [l^ which is located at Bq ~ 589 G and has a width 
of AB ~ 1.5 G [nl- 



The details of our experimental setup and procedure 



have already been described elsewhere [17|, |20|| and shall 
only be summarized here. We produce a BEC of ap- 
proximately 25,000 atoms about 10 G above the reso- 
nance where the scattering length is still close to its 
background value. Once the BEC is obtained, we adi- 
abatically shape the trapping potential to the desired as- 
pect ratio A within 25 ms. In order to be able to vary 
A over a large range, we generate the trapping poten- 
tial by a crossed optical dipole trap and a superimposed 
one-dimensional optical lattice along the z-direction (see 



Fig. IC). The two lattice beams (wavelength Aiatt = 
1064 nm, waist wiatt = 110 /xm, maximum power per 
beam Piatt = 5 W) propagate in the xz-plane under a 
small angle of ??/2 = 4° with respect to the x-axis. This 
configuration creates a standing wave along the z-axis 
with a spacing d = Aiatt/[2 sin(?9/2)] — 7.6 ^m. Due to 
the large spacing we load at most two sites when ramp- 
ing up the optical lattice. Tunneling processes are com- 
pletely negligible on the timescale of our experiments. By 
varying the powers in the beams we are able to provide 
nearly cylindrically symmetric traps, with aspect ratios A 
between ~ 1/10 and ^10, while keeping the avera ge t rap 
frequency Cu — {uj^uj^Y^^ approximately constant 2l| . 

We then ramp the magnetic field within 10 ms to adjust 
the value of the scattering length. The current provid- 
ing the magnetic field is actively stabilized on the 10~^ 
level [23 | , which results in a resolution of Aa ~ ag around 
the zero crossing of the scattering length. After an addi- 
tional holding time of 2 ms we finally switch off the trap 
and take an absorption image along the x-axis, after a 
time of fiight of 5 ms. The BEC atom number and radii 
are obtained by fitting the density profile using a bimodal 
distribution [20l| . Knowing the atom number and radii we 
can calibrate the scattering length a as a function of the 
magnetic field B [i3,[23|. 



We observe two effects when approaching the zero- 
crossing of the scattering length: The BEC shrinks in 
both directions due to the decreasing scattering length 
and the ellipticity of the cloud changes as a manifesta- 
tion of the enhanced dipolar effects [l7|. Finally, when 
we decrease the scattering length even further, the BEC 
atom number abruptly decreases. At this point the den- 
sity distribution does not show a bimodal shape any more 
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TABLE L Trap frequencies and aspect ratios of the traps 
that we used. The trap frequencies were measured by either 
exiting the center of mass motion or parametric heating and 
are accurate to about 10%. 



Trap 


i.^r/(27r) (Hz) 


uj.,I{2-k) (Hz) 


u)I{2-k) (Hz) 
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620 


0.11 
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400 
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330 
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720 
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but becomes thermal-like. The total atom number stays 
roughly constant during this collapse, excluding three- 
body loss processes causing the decrease in BEC atom 
number. The critical scattering length Ocrit where the 
condensate collapses depends strongly on the trap aspect 
ratio A (Fig. 2). For an isotropic trap (red) the collapse 
occurs at a ~ 15ao, whereas the pancake-shaped trap 
(blue) can even stabilize a purely dipolar BEC (a ~ 0). 
We repeated this experiment for all the six traps listed 



in Table m thereby covering a range of two orders of mag- 
nitude in the trap aspect ratio A. By fitting to the ob- 
served BEC atom numbers (Fig. 2) the threshold function 

iV = max[0,7Vo(a-acrit)''] , (5) 

where Nq^ amt and /3 are fitting parameters, we find 
the critical scattering length Ocrit • The simple functional 
form (Eq.[5]) was chosen because it accounts for the slowly 
decreasing BEC atom number when approaching the col- 
lapse point. The exponent (3 describing the steepness of 
the collapse was found to be /? ~ 0.2 for all traps. The ob- 
tained values of Ocrit versus the trap aspect ratio are plot- 
ted in Fig. 3 A. We observe a clear shift towards smaller 
a as A increases. For the most oblate trap (A = 10) we 
can reduce the scattering length to zero and hence access 
the purely dipolar regime experimentally. 

To get a more quantitative insight into the collapse 
threshold acrit(A) we numerically determine the critical 
scattering length (green curve in Fig. 3 A) as follows. 
The ground state wave function $(r) of a BEC can be 
found by the minimization of the Gross-Pitaevskii energy 
functional [l[ 



Em 



|^|V$p + Ftrapl^P + f 1*1^ + \\^? J t/dd(r - r')|$(r')pdr' 



dr. 



(6) 



Here the first term corresponds to the kinetic energy i?kin, 
the second to the potential energy i?pot in the trap, while 
the third represents the contact interaction energy i?cont j 
where g = A-Kh'^a/m is the coupling constant. The last, 
non-local term Edd arises from the magnetic DDI [3], 
where 



Udd{r) 



3 cos2 e 



4tt 



(7) 



is the interaction energy of two magnetic dipoles fi 
aligned by an external field. Here r is the relative po- 
sition of the dipoles and 6 the angle between r and the 
direction z of polarization. In order to obtain an estimate 
of flcrit we calculate the energy E(crr, Oz) of a cylindrically 
symmetric Gaussian wave function 

*(^' ^) = 1 ^3/2^9 _ ) exp 



1 



Ho 



with ar and Cz as variational parameters 
ansatz, where Oh 



total energy are the zero point fluctuations 



(8) 

Using this 
the contributions to the 



kin 



(9) 



the potential energy 



E, 



pot 



1 



2crf + X'cr 



NUco 4A2/3 
and the mean-field interaction energy 



Er. 



+ Edd Nadd 



1 



Nhui 



- 



(10) 



(11) 



Here /(k) is a monotonically decreasing function of the 
condensate aspect ratio k = cTr/cz with the asymptotic 
values /(O) = 1 and /(oo) = —2, arising from the non- 
local term in Eq. [S] [24|. The function / vanishes for 
K ~ \ implying that for an isotropic density distribution 
the magnetic DDI does not contribute to the total en- 
ergy. To obtain Ocrit we lower the scattering length until 
the energy landscape E{ar,<yz) does not contain a min- 
imum for finite ct^ and cTz any more (see Fig. 3 B-E). 
Starting with large values a > add we find that E{ar, (Jz) 
supports a global minimum for finite Ur and Uz indepen- 
dently of A and thus the BEC is stable (Fig 3 B). Going 
below a ~ add the absolute ground state is a collapsed in- 
finitely thin cigar-shaped BEC {ar — > 0) and the possible 
existence of an additional local minimum (corresponding 
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FIG. 3: (A) Stability diagram of a dipolar BEG. Experimental (green squares) and theoretical (green line) values of the critical 
scattering length ttcrit are plotted as a function of the trap aspect ratio. The theory curve is obtained for 20,000 atoms and 
an average trap frequency = 27r x 700 Hz (the average values we find for our six traps). The red curve (magnified in the 
inset) marks the stability threshold for a BEG with pure contact interactions using the same parameters. In grey we plot the 
asymptotic stability boundary (A'' — > oo) which for A — > (A ^ oo) converges to add (— 2add), see text. (B-E) Behavior of the 
energy landscape E{ar, (Jz). Lines of equal energy are plotted for fixed A = 10 and four different values of the scattering length 
a (blue dots in (A)). For Ocrit < a < add (C) the collapsed prolate ground state emerges [or ^ at finite a^) and the BEG 
becomes metastable. 



to a metastable state) is determined by the trap aspect 
ratio A (see Fig 3 C, where add > a > acrit and Fig D, 
where a = acrit)- Finally, below a ^ — 2add (Fig 3 E) 
the metastable state vanishes for any A and the BEC is 
always unstable [l^, [l3l ■ 

Considering the limit iVadd/aho ^ 1 where the terms 
(O and PH)) can be neglected [23| (grey curve in Fig. 3 A) 
we gain further insight into the nature of the dipolar col- 
lapse. In this case the stability is governed by the compe- 
tition between the contact and dipole-dipole interaction 
only, that is by the sign of the last term in Eq. [11] Hence 
the critical scattering length is implicitly given by 

acrit(A) = add/(K(A)). (12) 

The asymptotic behavior of the theory curve Ccrit = Odd 
for A ^ (respectively acrit = — 2add for A oo) now 
becomes apparent as for extremely prolate (respectively 
oblate) traps the cloud shape follows the trap geometry 
and / takes on its asymptotic values. Another particular 
point is acrit = marking the aspect ratio A needed to 
stabilize a purely dipolar BEC. More precisely, as /(I) = 
0, we search for the trap in which the ground state of 
a purely dipolar BEC is isotropic. As the DDI tends 
to elongate the BEC along the z-direction and shrink it 
radially, it is clear that the desired trap is oblate. Using 
our model we obtain the criterion A > Ac ~ 5.2 for a 
purely dipolar BEC to be stable, a result that agrees 
well with the values found in 

The grey curve in Fig. 3 A that we obtain by numer- 
ically solving Eq. [T^] shows a universal behavior in the 
sense that in the large N limit acrit(A) does not depend 



any more on the absolute values of the trap frequencies 
and N. This fact clearly distinguishes the dipolar col- 
lapse from the pure contact case (red curve in Fig. 3 A). 
The former is ruled by the interplay between contact 
and dipolar interaction whereas in the latter the zero 
point energy and the contact interaction rival against 
each other. Due to the different iV-scaling of the two com- 
peting terms in the pure contact case, the A-dependance, 
which is already weak for finite N [27|, completely van- 
ishes in the limit of large N as the stability criterion 
reads acrit(A) = (see Eq. [T]and red curve in Fig. 3 A). 
Furthermore the stability threshold obtained here applies 
for any dipolar system like e.g. hetcro-nuclear molecules, 
where the only difference is the specific value of add- 

In spite of the simplicity of our model we find good 
agreement between experiment and theory (Fig. 3 A). 
We checked that the different atom numbers and mean 
trap frequencies that we find for the six traps modify the 
green curve by much less than the error bars which arise 
mainly from the calibration of the scattering length. For 
the most oblate trap (A — 10) the 1/e- lifetime of the 
purely dipolar BEC (a — 0) decreases to ~ 13 ms. 

In conclusion, we experimentally mapped the stabil- 
ity diagram of a dipolar BEC. The dependance on scat- 
tering length and trap aspect ratio agrees well with a 
simple model based on the minimization of the energy 
of a Gaussian ansatz. By using a pancake-shaped trap 
we were able to enter the regime of purely dipolar quan- 
tum gases. This work opens up the route to new and 
exciting physics [2^. A clear subject for future studies is 
the dynamics of the dipolar collapse, which might show 
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anisotropic features. Another remarkable property of a 
dipolar BEC in a pancake-shaped trap is the existence of 
a roton minimum in its Bogoliubov spectrum [isj . Fur- 
thermore, close to the collapse threshold, the existence 
of structured ground states is predicted 28, |2§|, a pre- 
cursor for the supersolid phase [s^] that is expected to 
appear in dipolar BECs in three dimensional optical lat- 
tices. Finally, a field that has gained increasing interest 
in the recent past is the study of unusual vortex lattice 
patterns in rotating dipolar BECs [sij . 
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